† Corresponding author. E-mail:
We present the application of differential quadrature (DQ) method for the buckling analysis of nanobeams with exponentially varying stiffness based on four different beam theories of Euler–Bernoulli, Timoshenko, Reddy, and Levison. The formulation is based on the nonlocal elasticity theory of Eringen. New results are presented for the guided and simply supported guided boundary conditions. Numerical results are obtained to investigate the effects of the nonlocal parameter, length-to-height ratio, boundary condition, and nonuniform parameter on the critical buckling load parameter. It is observed that the critical buckling load decreases with increase in the nonlocal parameter while the critical buckling load parameter increases with increase in the length-to-height ratio.
Proper knowledge of mechanical properties is quite important in the design of various kinds of materials.[1–3] Due to their excellent physical, mechanical, and electrical properties,[4] nanostructures have attracted much attention among the scientists/researchers to develop innovatory applications in the field of nanomechanics. Proper understanding of their mechanical behavior is a key factor in the production of such engineering structures. Among these nanostructures, nanobeams attract more attention due to their great potential in engineering applications such as nanowires, nanoprobes, atomic force microscope (AFM), nanotube resonators,[5] nanoactuators,[6] and nanosensors. When nanostructure elements are subjected to compressive in-plane loads, these structures may buckle. Previously researchers have proposed that the surface stress itself can induce the buckling of a nanobeam with large deformation.[7] Figure
For a proper design of nanobeams, it is important to accurately predict the buckling characteristics of the nanobeams. The most challenging area in the field of nanomechanics is the study of structures at very small length scales. When the structures are of nanoscale size, classical continuum theories cannot accurately predict the mechanical behavior of such structures. Hence various size-dependent continuum theories such as strain gradient, couple stress, micropolar, and nonlocal elasticity have been developed to capture the size effects. Among these theories, the nonlocal elasticity theory proposed by Eringen[8] has gained much attention among the researchers because of its efficiency and simplicity to analyze the behavior of nanostructures. Applications of nonlocal continuum mechanics can be found in the areas of lattice dispersion of elastic waves, fracture mechanics, dislocation mechanics, wave propagation in composites, etc. It has been observed that the results obtained by the nonlocal continuum models are different from those previously obtained by the classical continuum theories, which indicates the influence of the size effects on the behavior of structures at nanoscale. Therefore, inclusion of the size effects in the buckling analysis of nanobeams is necessary. The presence of a nonlocal parameter in the constitutive relation indicates the inclusion of small scale effect. The nonlocal parameter has been calibrated using molecular dynamics to find its actual values.[9,10] The effect of the nonlocal parameter on the frequency and mode shapes is significant in the case of nanobeams.[11] Various relevant studies in the context of buckling analysis of nanobeams based on the nonlocal elasticity theory are shown in the subsequent paragraphs.
Buckling analysis of nano rods or tubes has been investigated analytically by Wang et al.[12] based on nonlocal Euler–Bernoulli and Timoshenko beam theories. Analytical results have been shown for four types of boundary conditions (simply-supported (SS), clamped-simply-supported (CS), clamped-clamped (CC), and cantilever (CF)). Reddy[13] reformulated various nonlocal beam theories such as Euler–Bernoulli beam theory (EBT), Timoshenko beam theory (TBT), Reddy beam theory (RBT), and Levison beam theory (LBT) based on the nonlocal elasticity theory and analytical results have been shown for SS nanobeams. Static nonlinear postbuckling response of nanobeams has been studied analytically by Emam[14] and the results show that as the nonlocal parameter increases, the critical buckling load reduces while the amplitude of buckling increases. Analytical methods have also been used in some of the studies[15–20] related to nanobeams. One of the numerical methods, the Rayleigh–Ritz method, has been employed in most of the studies such as the buckling analysis of single-walled carbon nanotubes in thermal environments,[21] bending, buckling, and vibration problems of nonlocal Euler beams,[22] buckling analysis of single-walled carbon nanotubes subjected to different boundary conditions,[23] free vibration of nonhomogeneous Timoshenko nanobeams,[24] free vibration of rectangular nanoplates,[25] etc. A differential transformation method has been applied by Pradhan and Reddy[26] to predict the buckling behavior of a single-walled carbon nanotube (SWCNT) embedded on a winkler foundation. Sahmani and Ansari[27] developed a state-space method to study buckling of nanobeams based on Euler–Bernoulli, Timoshenko, and Levison beam theories. The effect of temperature on the buckling analysis of SWCNTs embedded in an elastic medium has been investigated by Narendar and Gopalakrishnan.[28] The nonlocal shell model has been used by Yan et al.[29] to study the nonlocal effect on axially compressed buckling of triple-walled carbon nanotubes under the influence of temperature. Some of the studies in context of buckling of single-layered graphene sheets based on plate theories were presented in Refs. [30]–[33].
Structural members with variable cross section are frequently used in civil, mechanical, and aeronautical engineering to satisfy architectural requirements. In practical cases such as space structures, this type of buckling analysis plays an important role in design. Many engineers currently design light slender members with variable cross sections to construct ever-stronger and ever-lighter structures. Unfortunately, design engineers are lacking proper knowledge on the design of nonuniform structural elements since most of the design specifications are available for uniform elements. Hence there is a need for buckling analysis of nonuniform structural elements. The structural design must ensure the stability of the initial equilibrium shape to evaluate correctly the critical loads corresponding to neutral and unstable equilibriums. Many structural elements have variable flexural rigidity which may be due to different reasons (technological, economical, etc.) and the present paper studies buckling analysis of such nonhomogeneous nanobeams.
Since conducting experiments with nanoscale size is quite difficult, so efforts are always being made by the researchers to develop efficient numerical or analytical methods for obtaining better results. The literature reveals that analytical methods have been used in most of the studies. However, it is not always possible to obtain analytical solutions for complicated geometries. Researchers have developed a few numerical methods in the buckling analysis of nanobeams. One must have sound knowledge about variational principles for applying numerical methods such as finite element method and Rayleigh–Ritz method. Again, subsequent application of variational principles often requires a proper understanding of principles of mechanics. This has motivated the search for an approximate computational technique for the direct solution of problems without recourse to variational principles. In this context, one may use the differential quadrature method.
The differential quadrature (DQ) method was introduced by Bellman and Casti[34] for the first time. This is an efficient numerical method for the solution of linear and nonlinear partial differential equations. Then Bert et al.[35] applied the DQ method in structural problems. Since then, researchers have been investigating linear and nonlinear structural problems using the DQ method. Different procedures have been used by the authors to implement the boundary conditions in the DQ method. Firstly, the δ technique was proposed by Bert et al.[36] to implement the boundary conditions. In this procedure,[37] one boundary condition is implemented at the boundary point and the other boundary condition is implemented at a distance δ from the boundary point. There are two major drawbacks to this approach. Firstly, since the implementation of the boundary condition at the δ point is an approximation of the true boundary condition which should be implemented at the boundary, accurate numerical solutions lie on the smaller value of δ. Secondly, smaller value of δ causes the solutions to oscillate since the weighting coefficient matrices become highly ill-conditioned. The above mentioned approach is suitable for clamped ends but not suitable for simply-supported and simply-supported-clamped. Seeing demerits of this approach, Bert[38–42] proposed a new approach in applying the boundary conditions. In this approach, only one boundary condition is numerically implemented and the other boundary condition is built into the DQ weighting coefficient matrices. The following are some of the advantages of this approach. (i) The boundary conditions are properly satisfied since they are applied at the boundary points, therefore, better results may be obtained. (ii) The effect of δ on the results is eliminated. (iii) Excellent results are obtained with less computational effort.
On the other hand, the previous studies present investigations on the variable cross section in case of beams[43,44] and a few in case of nanostructures.[45] In case of buckling analysis of beams, the previous authors have used various ways of nonuniform cross section. In practical situations, stiffness of nanobeams may vary exponentially. To the best of the authors’ knowledge, till now, no work has been reported on buckling analysis of nanobeams having exponentially varying stiffness. It is also found that no literature on application of DQ in the buckling analysis of nanobeams subjected to different boundary conditions is available. In this article, the DQ method is employed to investigate nonuniform cross section on buckling analysis of nanobeams. Different beam theories such as EBT, TBT, RBT, and LBT are taken into consideration. Governing equations are converted into a single unknown function. Critical buckling loads are shown for various types of edge conditions such as SS,CS, CC, guided, and SS-guided boundary conditions. Effects of nonlocal parameter, nonuniformity, beam theories, boundary conditions, and length–height ratio have been investigated in this article.
This study is carried out on the basis of four beam theories: Euler–Bernoulli beam theory, Timoshenko beam theory, Reddy beam theory, and Levinson beam theory. The present formulation is based on the nonlocal elasticity theory of Eringen.[8] Displacement fields, nonlocal constitutive equations, and stress–strain relations of the beam theories are well known and may be found in Ref. [13]. The governing equations which are converted to a single variable form are presented here. One should note that x, y, and z coordinates are taken along the length, width, and thickness of the beams, respectively (Fig.
The governing equation in terms of displacement may be written as[13]
The governing equations for TBT nanobeams are given by[13]
In this case, the governing equations in terms of displacements may be written as[13]
By eliminating φ from Eqs. (
The governing differential equations may be written as[13]
Next, we assume an exponential variation of the flexural stiffness (EI) since the flexural stiffness of the SWCNT may not be constant for a geometrically non-uniform beam model of SWCNT. For this, we consider a beam with a nonuniform variation of the cross section along the length. Based on the proposed exponential variation, the flexural stiffness is defined as follows:
In this analysis, we have taken the following non-dimensional terms:
Below we show the non-dimentionalized forms of the governing differential equations for EBT, TBT, RBT, and LBT respectively:
The derivatives of displacement function
Let us now denote
For the simply supported boundary condition, equation (
Now, since
Similarly, for the fourth order derivative,
For the clamped-simply supported boundary condition, proceeding in a similar fashion as before, we have
For the clamped-clamped boundary condition, we have
For the guided boundary condition, we have
For the simply supported-guided boundary condition, we have
It is noted that when substituting the derivatives in the governing differential equations, one has to use
In this section, numerical results have been obtained by solving Eq. (
First, to validate the present method, we compare our results of the critical buckling load parameter with those available in the literature. For the validation, we consider uniform (
To find the minimum number of grid points for obtaining accurate results, a convergence study is carried out for EBT and TBT nanobeams. To show how the solution is affected by the grid points, the variation of the critical buckling load parameter with the number of terms (N) is shown in Fig.
In this subsection, the significance of the scale coefficient is highlighted. First we define the buckling load ratio as
In this subsection, the effect of non-uniformity η on the critical buckling load parameter is illustrated. It may be noted that η takes positive values only and different value of η gives different cross section of the nanobeams. This analysis will help design engineers to have an idea of the values of the critical buckling load parameter. Figure
Another important factor that design engineers should keep in mind is the length-to-height ratio. To investigate the effect of length-to-height ratio
Modeling of nanostructures based on beam theories is one of the important areas in the field of nanotechnology. Various beam theories used in this investigation are discussed below. The first and simplest beam theory is EBT. In this theory, the transverse shear and transverse normal strains are ignored. The next theory developed in the hierarchy of beam theories is TBT. According to this theory, a constant state of transverse shear strain with respect to the thickness coordinate is assumed. Therefore, a shear correction factor is required in order to compensate the error due to the above assumption. In RBT and LBT, the transverse shear strain varies quadratically and also vanishes on the top and bottom planes of the beams. Thus the shear correction factor is not taken into consideration. EBT and TBT are the first order beam theories, while RBT and LBT are the third order beam theories. To investigate the effect of various beam theories such as EBT, TBT, RBT, and LBT on the buckling load parameter, the variation of the critical buckling load parameter with scale coefficient is shown in Fig.
For designing engineering structures, one must have proper knowledge about the boundary conditions. It will help engineers to have an idea without carrying out detailed investigation. Therefore, analysis of the boundary conditions is quite important. In this subsection, we consider the effect of the boundary condition on the critical buckling load parameter. Figure
Buckling analysis of nonuniform nanobeams based on four different beam theories including EBT, TBT, RBT, and LBT has been carried out. The nonuniform cross section is assumed by taking exponentially varying stiffness. The differential quadrature method has been employed and the boundary conditions have been substituted in the coefficient matrix. New results have been shown for two types of boundary conditions, i.e., guided and SS-guided boundary conditions. The numerical results are presented to show the effects of the nonuniform parameter, the nonlocal parameter,